全模颤振风洞试验三索悬挂系统多体动力学分析

全模颤振风洞试验需要通过软支撑系统模拟飞行器的自由飞行状态并调整模型姿态达到配平状态。参考NASA双索悬挂方案，提出了一种两电机驱动的三索悬挂系统，利用后方两索的同向/反向联动实现模型俯仰和滚转姿态的调整，利用弹簧刚度以及钢绳张力设计实现支撑频率要求。基于柔性多体动力学方法，建立了包括飞行器刚体模型、柔性索、滑轮、弹簧、气动力模型、伺服电机控制在内的复杂系统动力学模型，其中，利用任意拉格朗日-欧拉（ALE）变长度索单元描述钢绳，利用不约束物质坐标的索结点约束描述钢绳与滑轮相互作用，利用索结点物质输运速度约束描述伺服电机绞盘，利用飞行力学的气动力模型描述吹风下的气动力。基于该模型，通过小扰动响应辨识研究了弹簧刚度、钢绳张力、连接点位置等因素对支撑频率的影响规律，并分析了系统姿态调整能力，俯仰调整范围达到-12.5°~12.5°，滚转调整范围达到-45°~45°。采用滑轮处电位计测量的钢绳相对位移作为反馈信号，基于设计的控制律利用多体动力学求解器与Simulink对风洞吹风下的姿态调整过程进行仿真，模型达到配平状态，获得了吹风下的索拉力和伺服电机功率，为系统设计提供基础。     

前缘缝翼构型平尾直升机气动特性分析

使用数值模拟方法研究了前缘缝翼翼型气动特性,并将该前缘缝翼翼型应用于直升机平尾,分析了该直升机全机气动特性,最后通过风洞试验对数值模拟结果进行验证。结果表明,前缘缝翼构型平尾能有效改善直升机的纵向静稳定性。     

全模颤振双索悬挂系统刚体模态频率研究

研究绳系结构和绳索预紧力对全模颤振双索悬挂系统刚体模态频率的影响。首先，基于绳索并联机构理论，对双索悬挂系统进行静力学建模，引入加权矩阵，推导了双索悬挂系统静刚度模型；其次，建立双索悬挂系统无阻尼振荡方程，分析双索悬挂系统刚体模态频率变化规律；进而，搭建了双索悬挂系统地面样机，开展系统模态频率测试试验，研究系统刚体模态频率的影响因素及其变化规律；最后，参考试验结果修正所建数学模型。结果表明，随着绳索预紧力的增加，系统刚体模态频率呈上升趋势，但各阶刚体模态的固有频率上升速率不同，平动模态的固有频率受绳索预紧力影响比转动模态小得多；在转动模态方面，滚转模态的固有频率最大且上升快；在平动模态方面，升沉模态的固有频率大于横侧滑模态的固有频率；双索经过的机身前后滑轮与飞机模型质心之间的相对位置变化会不同程度地影响支撑系统俯仰模态、偏航模态的固有频率，但几乎不影响其滚转模态的固有频率。研究发现，通过控制悬挂绳索预紧力的大小，合理设计机身上的滑轮安装位置，能够有效降低双索悬挂系统刚体模态频率。研究结果对于优化设计全模颤振双索悬挂系统具有指导意义。     

翼型低雷诺数层流分离现象随雷诺数的演化特征

层流分离现象是翼型低雷诺数条件下出现的典型流场特征。层流分离流动中包含流动分离、转捩、再附等非定常流动结构，层流分离流动的形成与演化会对翼型气动特性产生恶化作用。采用大涡模拟（LES）方法对低雷诺数范围内不同雷诺数下的翼型层流分离流动开展精细数值模拟，研究了雷诺数对翼型气动特性的影响规律及作用机理。LES方法采用隐式亚格子模型，基于结构化拼接网格，对流项离散和时间推进方法分别采用AUSM⁺格式以及双时间步方法。验证算例计算结果表明数值模拟方法的正确性及可靠性，雷诺数对翼型气动特性具有显著影响。随雷诺数降低，时均分离泡外形增大、位置后移，平均阻力系数增大，特别是在较低雷诺数下，翼型升阻力系数随时间出现振荡现象。进一步研究表明，造成不同时均分离泡形态和气动特性的原因在于翼型上表面分离剪切层的失稳与转捩特征。随雷诺数降低，流动黏性增大，导致分离剪切层速度梯度减小，流动发生转捩及再附位置后移，直至翼型表面不再发生转捩和再附。     

不同后掠角大展弦比复合材料机翼气动特性

为了研究不同后掠角复合材料机翼的气动特性,建立了载荷模型;采用 K方法通过 Theodorsen函数推导得到了机翼颤振速度的计算公式,分析了 0°、5°、10°、15°、20°、25°和 30°共 7种不同后掠角机翼的模态,得到了模态频率和模态振型;仿真计算得到了机翼的颤振速度和发生颤振时翼尖的最大垂直位移。结果显示:机翼后掠角为 15°时,其颤振速度和翼尖最大垂直位移均在安全范围内,因此,后掠角 15°,展弦比为 26的机翼为飞行器机翼设计最优值。     

Stall flutter prediction based on multi-layer GRU neural network

The modeling of dynamic stall aerodynamics is essential to stall flutter, due to the flow separation in a large-amplitude pitching oscillation process. A newly neural network based Reduced Order Model (ROM) framework for predicting the aerodynamic forces of an airfoil undergoing large-amplitude pitching oscillation at various velocities is presented in this work. First, the dynamic stall aerodynamics is calculated by solving RANS equations and the transitional SST-γ model. Afterwards, the stall flutter bifurcation behavior is calculated by the above CFD solver coupled with structural dynamic equation. The critical flutter speed and limit-cycle oscillation amplitudes are consistent with those obtained by experiments. A newly multi-layer Gated Recurrent Unit (GRU) neural network based ROM is constructed to accelerate the calculation of aerodynamic forces. The training and validation process are carried out upon the unsteady aerodynamic data obtained by the proposed CFD method. The well-trained ROM is then coupled with the structure equation at a specific velocity, the Limit-Cycle Oscillation (LCO) of stall flutter under this flow condition is predicted precisely and more quickly. In order to predict both the critical flutter velocity and LCO amplitudes after bifurcation at different velocities, a new ROM with GRU neural network considering the variation of flow velocities is developed. The stall flutter results predicted by ROM agree well with the CFD ones at different velocities. Finally, a brief sensitivity analysis of two structural parameters of ROM is carried out. It infers the potential of the presented modeling method to depict the nonlinearity of dynamic stall and stall flutter phenomenon.     

重型燃气轮机高雷诺数CDA叶型转捩特性数值计算

为研究重型燃气轮机的压气机叶片在高雷诺数工况下的气动性能,基于Gamma-Theta转捩模型的雷诺时均方程对某可 控扩散叶型进行了数值计算。通过对比不控制马赫数与控制马赫数,分析高雷诺数对可控扩散叶型气动性能及转捩特性的影响。 结果表明：在不控制马赫数条件下,在零攻角时,雷诺数从7×10 5 增大为9×10 5 ,总压损失增加了约391.95%;在高雷诺数工况下随 着雷诺数的增大,叶片流动损失不断增大,叶片可用攻角范围减小,同时在叶片吸力面出现激波,干扰转捩的产生。在控制马赫数 条件下,当Ma=0.6时,在零攻角工况下,雷诺数从8.2×10 5 增大为1×10 7 ,总压损失减小了约38.98%,吸力面转捩起始点从4.78%弦 长处前移至1.11%弦长处;在高雷诺数工况下,叶片流动损失随着雷诺数的增大不断减小,吸力面转捩位置前移。     

低雷诺数仿生分离流翼型气动特性初探

为了探索适合低雷诺数微型飞行器的翼型形式,基于对自然界鸟类和昆虫滑翔飞行时翅膀形状的观察,设计出一种由前缘削尖平板和后缘圆弧翼型组合而成的仿生分离流翼型。数值研究结果表明,气流在削尖平板的前缘点强制分离,形成大范围低压分离流动,随后在后部圆弧翼上表面再附形成稳定低压涡流区,从而实现较高的气动效率和较强的抵抗大气湍流的能力。上削尖平板可以使流动分离点固定在削尖点。相对于单独平板,仿生分离流翼型的升力系数有大幅提高,迎角为4°时提高了112%。此外,仿生分离流翼型可以在较宽的迎角范围内（4°~20°）保持高升力,但是迎角增加,阻力也快速增大,因此小迎角情况下（小于4°）气动效率更优。      

适用于多迎角结冰试验的混合翼型设计

为突破结冰风洞对翼型模型尺寸的限制,提出了一种新的混合翼型设计方法,可使用一套混合翼型模拟原始翼型在不同迎角下的结冰试验,弥补了以往混合翼型只能在单个设计迎角下使用的缺陷。方法采用多段翼的形式设计混合翼型,以多目标迎角等结冰试验条件作为设计输入,优化设计主翼外形和襟翼的位置、偏转角度,利用襟翼位置和偏转角度的变化实现混合翼型和原始翼型前缘表面的压力系数在不同迎角下能够保持一致,继而保证前缘结冰外形的一致性。设计的混合翼型较原始翼型弦长减小接近40%,在冰风洞中对混合翼型和原始翼型在目标结冰条件下进行试验,对比结果显示,在选取的多个目标迎角下,混合翼型和原始翼型二维截面的结冰外形基本一致,提出的混合翼型设计方法能够有效的模拟原始翼型在不同迎角下的结冰外形。     

惯性参数对飞翼体自由度颤振特性的影响规律研究

飞翼飞行器刚体短周期模态频率较高,易与一阶弹性弯曲模态耦合发生一种特殊的颤振——体自由度颤振。采用风洞实验与频域计算相结合的手段,开展了惯性参数(俯仰转动惯量、质心位置)对飞翼飞行器体自由度颤振特性(颤振速度和颤振频率)的影响规律研究。实验和计算结果表明:俯仰惯量和质心位置会明显改变体自由度颤振频率与速度,颤振实验与计算结果一致性较好。俯仰惯量增加,颤振频率降低,沉浮约束时颤振速度基本不变,沉浮自由时颤振速度增大;质心位置前移,俯仰模态频率与阻尼同时增加,俯仰与一阶弹性弯曲模态耦合更容易,但发散模态分支的阻尼也更大,导致颤振速度先降低后增大,颤振频率单调增加。